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ABSTRACT 


The  amplitude  and  power  of  a  large  family  of  radio  signals 
are  observed  to  have  log-normal  probability  density  func¬ 
tions.  Among  these  are  signals  propagated  through  random 
inhomogeneous  media,  a  notable  example  being  low  frequency 
atmospheric  radio  noise.  Of  greater  importance  are  certain 
radar  targets  that  have  been  observed  to  have  essentially 
log-normal  density  functions.  Both  ships  and  space  vehicles 
may  fall  into  this  category.  Curves  of  probability  of  detection 
versus  signal-to-noise  ratio  for  the  case  of  log-normal  signals 
in  Gaussian  noise  have  been  computed  and  are  presented  in 
this  paper.  The  curves  apply  for  square-law  detection  with 
varying  degrees  of  post-detection  linear  integration.  Both 
slowly  and  rapidly  fluctuating  signals  are  considered.  It  is 
shown  that  for  log-normal  signal  distributions  having  large 
variances  the  probability  of  detection  differs  significantly 
from  that  obtained  using  curves  based  on  an  assumed  Rayleigh 
signal  distribution. 
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I.  INTRODUCTION 


The  detection  of  fluctuating  signals  in  the  presence  of  Gaussian  noise 
has  been  studied  widely  [l]-[4]  since  the  extensive  pioneering  treatment  of 
detection  theory  for  constant  amplitude  pulsed  signals  in  noise  by  Marcum  [5j. 
All  of  these  studies  are  extensions  of  Marcum's  work.  In  almost  all  cases, 
the  fluctuating  signal  amplitude  has  been  assumed  to  be  Rayleigh  distributed. 
This  choice  is  due  to  the  well-known  fact  that  the  sum  of  a  large  number  of 
independent  random  vectors  having  uniform  phase  distributions  is  a  Rayleigh 
vector,  i.e. ,  one  having  a  uniform  phase  distribution  and  a  Rayleigh  ampli¬ 
tude  distribution,  provided  only  that  the  number  of  contributing  vectors  is 
very  large  and  that  no  single  vector  contributes  significantly  to  the  total  power 
[6].  Swerling  [1]  has  presented  curves  of  the  probability  of  detection  as  a 
function  of  the  signal-to-noise  ratio  for  both  Rayleigh-distributed  signals  and 
for  signals  whose  power  is  distributed  Chi-square  with  four  degrees  of  free¬ 
dom.  Since  N,  the  number  of  pulses  integrated,  is  presented  as  a  parameter, , 
and  since  a  square-law  detector  is  assumed,  results  for  signals  whose  power 
is  distributed  Chi-square  with  2N  degrees  of  freedom  may  be  inferred.  It  is 
only  necessary  to  observe  that  the  sum  of  a  Rayleigh  noise  vector  and  a 
Rayleigh  signal  vector  is,  itself,  a  Rayleigh  vector  and  that  the  post -detection 
sum  of  N  independent  signal  plus  noise  variates  is  then  the  sum  of  the  squares 
of  2N  independent  Gaussian  variates  and,  hence,  is  Chi-square  with  2N  degrees 
of  freedom.  As  Swerling  has  noted  [2],  one  can  construct  a  wide  variety  of 
probability  density  functions  from  the  Chi-square  family.  In  choosing  a 
member  of  the  Chi-square  family  to  fit  a  particular  distribution,  one  is  free 


k>  '  (wii  I  >.  i  r.um-lrrs,  e  g.  ,  th>  mean  anti  the  number  of  degrees  of  free- 

Utile  rent  degrees  ol  skewness  are  obtained  as  the  number  of  degrees 

•  >!  freedom  is  varied. 

i  lie  assumption  of  signals  composed  of  a  large  number  of  small  inde¬ 
pendent  random  contributions  appears  to  be  reasonable  in  many  cases,  e.g., 
Hie  signal;,  scattered  by  many  radar  targets  consist  of  contributions  from 
mimei  uu..  and  diverse  reflecting  elements.  It  has  become  common  practice 
io  treat  most  fluctuating  radar  returns  as  having  Rayleigh-distributed  ampli¬ 
tude:,.  The  experimental  evidence  for  the  existence  of  Rayleigh-distributed 
radar  returns  is  not  altogether  conclusive,  however.  Early  investigators  [7] 
lotind  that,  while  target  cross  section  distributions  followed  the  Rayleigh  law 
reasonably  well  for  low  cross  sections,  they  often  exhibited  a  pronounced 
high  cross  section  tail,  t.e.  ,  higher  probability  density  at  the  high  cross 
sections  than  indicated  by  the  Rayleigh  l.<w.  This  difference  is  consistent 
with  the  idea  that  a  few  large  surfaces  predominate  in  the  target.  The  case 
ol  ant  raft  appears  te  he  a  notable  exception  where  experimental  data  are 
i  oiu  lusively  Rayleigh.  There  hits  been  some  disagreement  on  this  point  [8j, 
but  Hie  major  body  of  opinion  supports  the  Rayleigh  interpretation  [9j,  [10]. 
More  recently,  Harrison  [11J  in  analysing  the  statistics  of  the  radar 
returns  from  ships  has  found  that  the  return  signal  strength  die? ribution 
i  losely  follows  a  log-normal  law.  Space*  vehicle  and  satellite  cross  sections 
have  also  typically  exhibited  non-Rayleigh  statistics.  They  are  usually 
characterised  by  a  very  large  dynamic  range  and  variance  [l^J-  Kennedy 
IS,  I.,;.-,  obtained  statistical  data  on  the  radar  cross  section  of  a  satellite, 

•  examining  the  distribution  over  aspect  angle  and  frequency.  The  resulting 
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cross  section  probability  density  functions  are  shown  in  Fig.  1.  One  of  the 
density  functions  was  obtained  by  measuring  over  aspect  angle  alone  with  all 
aspect  angles  assumed  to  be  equiprobable.  The  same  aspect  measurement 
procedure  was  followed  in  obtaining  the  other  density  function,  with  frequency 
averaging  being  obtained  by  stepping  the  frequency  of  alternate  radar  pulses. 
The  density  functions  appear  to  be  log-norrrial  and,  significantly,  exhibit 
mean  values  about  10  dB  above  the  median  values.  The  high  ratio  of  mean- 
to-median  values  is  due  to  the  pronounced  high  cross  section  tail  of  the 
distribution.  It  is  not  possible  to  obtain  such  a  mean- to- median  ratio  with 
any  member  of  the  Chi-square  family  suggested  by  Swerling. 

Weinstock  [)4]  has  studied  the  modeling  of  certain  classes  of  satellite 
bodies  and  has  attempted  to  fit  empirical  data  to  Chi-square  models.  The 
difficulty  in  fitting  members  of  the  Chi-square  family  is  evident  when  one 
considers  the  variance  and  mean-to- median  ratio.  Starting  with  Swerling1  s 
generalized  distribution  [2] 

=rro(f)k'lwtK'f)i 

we  have  for  the  variance 

-z 

_2.  .x 
D  (x)  = 

where  k  is  one-half  the  number  of  degrees  of  freedom  of  the  Chi-square 
distribution.  The  median,  x^,  is  given  by  the  solution  of 
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or 
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kx  /x 
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k-1 

u 


exp(-u)  du  = 


1 
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The  latter  expression  is  recognized  as  a  form  of  the  incomplete  gamma 

function.  Solving  for  x/x  with  k  =  1,  the  result  is  x/x  -  1.44.  For 

m  m 

k  =  2,  it  is  1.18,  and  it  approaches  unity  with  increasing  k.  Since  large 
variances  and  mean-to- median  ratios  may  characterize  radar  targets  of 
interest,  the  Chi-square  family  with  integral  k  may  be  unsuitable. 
Weinstock  [13  ]  has  allowed  k  <  1  in  order  to  fit  the  variance  and  mean 
to  empirical  data  but  has  found  disagreement  of  higher  order  moments 
and  a  lack  of  fit  on  the  tails  of  the  distributions.  He  found  that  for  a 
class  of  bodies  of  the  conducting  cylinder  type  a  reasonable  fit  could 
be  obtained  by  varying  k  between  0.  3  and  0. 7. 

Current  system  analysis  and  design  practice  uses  the  Swerling 
Chi-  squared  models  almost  exclusively,  with  the  mean  characterising 
the  distributions  treated.  Both  Weinstock  [14]  and  Swerling  have 
observed  that  if  a  single  parameter  is  to  be  used  to  characterize  a  distri¬ 
bution  the  median  is  preferred  over  the  mean.  The  median  is  less 
influenced  by  high  cross  section  tails  of  distributions  and  leads  to  a 
more  conservative  choice  of  design  parameters.  To  conclusively  check 
the  validity  of  this  assertion,  it  is  necessary  to  compute  the  probability 
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of  detection  as  a  function  of  median  signal-to- noise  ratio  for  distributions 
having  large  ratios  of  mean-to- median  values.  The  log-normal  distribution, 
in  addition  to  being  one  that  is  found  to  occur  in  practice,  is  one  that  repre¬ 
sents  an  extreme  in  high  cross  section  tails.  It  is  therefore  suitable  for 
checking  the  accuracy  of  design  procedures  that  use  existing  cross  section 
models  but  characterize  distributions  by  their  medians. 

This  paper  presents  the  results  of  a  numerical  computation  of  the 
probability  of  detection  for  a  large  family  of  fluctuating  signal  distributions 
that  are  log-normal. 

II.  LOG-NORMAL  SCATTERING  MECHANISMS 

The  nature  of  the  scattering  mechanisms  that  lead  to  log-normal 
density  functions  for  radar  targets  is  not  understood,  and  no  attempt  at  a 
complete  analysis  is  made  in  this  paper.  Nevertheless,  a  brief  conjecture 
on  this  point  seems  in  order. 

If  a  random  signal  amplitude  has  a  log-normal  distribution,  it  can 
be  represented  by 

x  =  A  e"  ^ 

where  t  i®  normally  distributed.  A  log-normal  distribution  occurs  in  radio 

propagation  through  inhomogeneous  media,  the  logarithmic  variable  t  being 

» 

interpreted  as  the  sum  of  independent  attenuation  factors  associated  with 
statistically  independent  regions  or  "blobs"  in  the  medium  [15].  In  the 
case  of  the  radar  cross  section  measurements  referred  to  above,  however. 
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it  seems  clear  that  the  signal  fluctuations  are  due  to  statistical  effects 
associated  with  the  targets.  One  possible  explanation  for  log-normal 
target  statistics  has  been  suggested  by  DuWaldt  [16]  who  noted  that  the 
effect  of  surface  roughness  on  the  specular  return  from  a  flat  plate  is 
to  reduce  the  cross  section  relative  to  that  of  a  smooth  plate  by  a  factor 
e  \  where  'll  is  proportional  to  the  mean- square  surface  depth  fluctuation. 

If  the  radar  return  at  any  aspect  angle  is  that  from  a  single  surface  selected 
from  an  ensemble  having  a  normally  distributed  degree  of  roughness,  log 
normal  amplitude  statistics  result.  Since  the  return  from  a  ship  at  any 
given  aspect  angle  is  largely  due  to  a  single  predominant  dihedral  corner, 
the  log-normal  character  of  ship  returns  might  thus  be  explained. 

f 

m.  THE  PROBABILITY  OF  DETECTION  FOR 
LOG-NORMALLY  DISTRIBUTED  SIGNALS 

The  computational  procedures  are  based  on  the  receiver  models  used 
by  Marcum  and  Swerling,  i.  e. ,  a  linear  receiver  preceding  a  square-law 
detector  and  followed  by  post-detection  integration  which  gives  equal  weight 
to  all  the  pulses  added.  As  in  Swerling  [1],  two  extremes  of  pulse-to-pulse 
correlation  are  treated.  These  are: 

Case  1:  All  pulses  occurring  during  an  integration  period  are 
perfectly  correlated.  Pulses  in  successive  integration 
periods  are  independent.  In  the  case  of  a  scanning  radar, 
this  implies  independent  fluctuations  from  scan  to  scan  but 
with  no  fluctuation  during  a  scan. 


Case  2:  The  fluctuations  are  independent  from  pulse  to  pulse. 
The  probability  density,  f^(y),  of  N  integrated  signal  plus  noise 
variates  at  the  detector  output  is  integrated  to  obtain  the  probability  of 
detection. 


’D(*)  =  f  fN(y’  dy 
Yb 


The  density  function  of  signal  plus  noise  and  the  probability  of  detection 
are  functions  of  some  measure  of  power  in  the  signal  distribution.  The 
average  signal-to-noise  power  ratio,  x,  is  used.  The  normalization  is  the 
same  as  that  used  by  Marcum,  and,  hence,  the  detection  threshold  level, 
Y^,  is  as  given  therein,  namely,  the  solution  of  * 


/uN/n  /*» 

1  -(l)  =/  fN<y’0)dy 

•'v 


where  n  is  the  number  of  pulses  in  the  time  required  for  the  false  alarm 
probability  to  reach  0. 5.  Values  of  Y^  as  a  function  of  N  and  n  are 
tabulated  in  Table  I, 

In  Case  1,  the  density  function,  f^(y,x),  is  obtained  by  averaging 
Marcum's  density  function  for  a  constant  signal  plus  noise  over  the 
distribution  of  the  signal-to- noise  power  ratio.  Thus, 


fN(y.x) s  f  *N(y/*)  <*x 


where  f^(y/x)  is  the  conditional  probability  density  of  the  integrated  detector 
output  signal  plus  noise  variate  y,  given  a  signal-to-noise  power  ratio  x  and 
where  f(x,x)  is  the  density  function  for  the  signal-to-noise  power  ratio.  For 
a  log-normally  distributed  signal, 

i_  expf-  (In  * 

•Jza  ox  L  2a 

2 

where  o  is  the  variance  of  In  x. 

For  the  presentation  of  the  results,  it  is  more  convenient  to  charac¬ 
terize  f(x,x)  by  the  parameters  x  and  p,  the  ratio  of  mean  to  median  values. 
In  terms  of  o  and  In  x, 

x  =  exp^g-  +  In  x  j 
p  =  exp(^-) 

In  case  2,  f^(y,  x)  is  the  result  of  an  N-fold  self- convolution  of 
f^y.x)  as  given  by  Eq.  (3). 

Details  of  the  numerical  computation  are  given  in  the  appendix. 

IV.  RESULTS  AND  CONCLUSIONS 

Curves  of  the  probability  of  detection  as  a  function  of  average  signal- 
to-noise  ratio  are  presented  in  Figs.  2-28  for  Cases  1  and  2,  with  the 
number  of  pulses  integrated,  N,  as  a  parameter.  Fixed  values  of  p  and  the 
false  alarm  number  n  were  used  in  plotting  the  curves,  namely. 


p  =  1,  2,  4,  8,  16  and  n  =  10^,  10®,  10^.  Figures  2-4  are  identical  to 
those  in  Marcum  since  for  p  =  1  the  variance  is  zero  and,  hence,  the  target 
is  not  fluctuating.  Cases  1  and  2  are  thus  identical  for  p  =  1,  In  Figs. 

29  -  43,  the  curves  have  been  replotted  with  p  as  a  parameter.  These 
figures  are  useful  for  interpolation  in  p. 

To  improve  the  accuracy  of  interpolation  in  p,  it  may  be  advantageous 

to  first  plot  the  curves  with  p  as  a  parameter  as  a  function  of  the  median 

signal-to-noise  ratio  (or  to  normalize  R  to  unity  median  signal-to-noise 
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ratio).  A  set  of  such  curves  for  n  =  10  is  presented  in  Figs.  44  -  50. 

The  curves  for  Case  1  intersect  near  50-percent  probability  of  detection 
and,  hence,  facilitate  interpolation.  A  very  different  behavior  is  observed 
in  Case  2  when  the  curves  are  plotted  versus  the  median  signal-to-noise 
ratio.  The  point  at  which  the  curves  intersect  shifts  to  higher  values  of 
probability  as  the  number  of  pulses  integrated  is  increased. 

The  form  of  presentation  is  that  of  Marcum  and  Swerling.  While  it  is 
recognized  that  this  format  is  not  always  convenient,  it  seems  preferable 
in  view  of  the  wide  circulation  of  the  curves  compiled  by  the  aforementioned 
authors  and  the  resulting  ease  of  comparison.  A  secondary  reason  is  the 
compression  of  the  curves  afforded  by  plotting  probability  of  detection 
versus  normalized  radar  range  which  is  inversely  proportional  to  the 
quarter-power  of  the  signal-to-noise  ratio.  This  compression  is  desirable 
because  of  the  flattening  of  the  curves  with  increasing  variance  of  the 
signal  distribution.  A  nonlinear  dB  scale  for  the  signal-to-noise  ratio  is 
added  for  convenience. 


As  p  is  increased  from  unity  at  constant  average  signal-to-noise  ratio, 
the  probability  of  detection  is  reduced  at  high  signal-to-noise  ratios  and 
increased  at  low  signal-to-noise  ratios.  The  sensitivity  to  a  change  in  p  is 
seen  to  be  greatest  for  small  p.  As  p  is  increased  to  very  large  values, 
the  probability  density  function  approaches  a  delta  function  at  the  origin. 

To  show  this,  it  is  first  noted  that  the  mode  of  the  log-normal  distribution 
is  at 


xmode 


The  peak  value  of  the  density  function  is 


Hence  • 


2/rrln  p  x 


!£  * 0 


It  can  be  shown  that  unit  weight  it  associated  with  the  delta  function  on 
examination  of  the  probability  that  the  signal-to-noise  ratio  exceeds  some 
positive  value  c 


10 


P(x>e )  =  f  *(x,x) 
•'e 


dx 


=  ±f 


»  2 
e  Z  dz 


where  z  =  (In  x  -  In  xJ/'/iZ  cr,  and 


E  -  *n  e  -  In  x  +  In  p 
V2  In  p 


Then,  if  e  is  fixed  at  any  finite  positive  value, 

lim  P(x>e)'=  0 

p-*« 

As  a  consequence,  the  probability  of  detection  for  any  finite  average  signal -to- 
noise  ratio  approaches  zero  as  p  is  increased  without  limit. 

The  relative  insensitivity  of  the  probability  of  detection  to  a  change  in 
p  at  p  =  16  is  due  to  the  fact  that  the  mode  of  the  distribution  is  already  36  dB 
below  the  average  value.  Hence,  the  probability  weighting  in  the  vicinity  of 
the  mode  is  not  significant  for  all  but  the  very  highest  signal-to-noise  ratios. 
At  larger  values  of  p.  the  probability  of  detection  is  governed  by  the  relatively 
flat  and  slowly  decreasing  tail  of  the  density  function.  The  behavior  of  f(x,  1) 
with  increasing  p  in  the  range  of  values  used  for  the  calculations  is  illustrated 
in  Figs.  61  and  62.  These  curves  are  indicative  of  the  behavior  for  all  x 
since 

1 


f(x,kx)  * 


f 


The  effect  of  variation  in  x  at  constant  p  is  illustrated  in  Fig.  63. 

It  is  of  interest  to  compare  the  results  with  those  obtained  by  Swerling 
using  Rayleigh-distributed  signal  amplitude.  Figures  51  -  56  show  the 

g 

comparison  for  n  =  10  and  p  =  2,  16  (Case  1)  and  p  =  2,  4,  8,  16  (Case  2). 
For  p  =  2,  there  is  relatively  good  agreement  between  the  curves  over  a  wide 
range  of  signal-to-noise  ratios.  (It  should  be  recalled  that  for  the  Rayleigh 
distribution  p  equals  1.44.  )  The  agreement  is  particularly  good  for  Case  2 
when  a  large  number  of  pulses  is  integrated.  This  is  to  be  expected  since 
the  sum  of  a  large  number  of  independent  variates  approaches  a  normal 
distribution  regardless  of  the  distribution  of  the  individual  variates.  The 

t 

curves  differ  widely  for  p  =  16,  however,  thereby  illustrating  the  in¬ 
applicability  of  Swerling's  curves  when  the  mean-to-median  signal  power 
ratio  is  very  large.  The  probability  of  detection  for  Casr  2  based  on  the 
normal  approximation  to  the  density  function  of  the  detector  output  is  also 
shown  in  Figs.  53  -  56.  The  details  of  the  computation  are  given  in  the 
appendix. 

Figures  57  -  60  show  the  above  data  plotted  versus  the  median  signal- 
to-noise  ratio.  For  the  curves  from  Swerling,  the  signal -to -noise,  ratio 
is  the  average  signal-to-noise  power  ratio.  Thus,  Rq  is  to  be  interpreted 
as  that  range  for  which  the  log-normal  target  has  unity  median  signal-to- 
noise  ratio  and  the  Rayleigh  target  has  unity  average  signal-to-noise  ratio. 
These  curves  permit  judgment  of  the  efficacy  of  entering  Swerling's  curves 
with  median  signal-to-noise  ratio  in  the  hope  of  obtaining  a  conservative 
estimate  of  performance.  It  is  seen  that  a  conservative  estimate  is  obtained 
for  p  »  2  in  both  Cases  1  and  2.  For  larger  values  of  p,  the  curves  for  the 


log-normal  distribution  fall  below  those  for  the  Rayleigh  distribution  over  a 
portion  of  their  range.  Nevertheless,  conservative  estimates  are  obtained 
almost  everywhere  for  p  =  16  in  Case  2.  In  Case  1  for  p=  16,  the  blip-scan 
ratio  predicted  by  the  Rayleigh  curves  never  exceeds  the  actual  one  by  as 
much  as  10  percent. 

When  the  statistical  description  of  a  radar  target  is  obtained  by 
averaging  over  all  of  the  aspect  angles  and  assuming  them  to  be  equally 
probable,  some  discretion  must  be  applied  in  using  the  resulting  model. 

It  is  entirely  possible  that  a  high  cross  section  tail  in  the  density  function 
results  from  specular  returns  at  relatively  few  aspect  angles.  If,  then,  the 

target  attitude  is  relatively  stable  during  a  radar  intercept,  the  equally 

'*1 

probabfe  aspect  assumption  can  lead  to  very  erroneous  results.  Conserva- 

* 

tively,  it  may  be  best  to  neglect  the  high  cross  section  tail  of  the  density 
function  when  a  lower  bound  on  detection  probability  is  desired.  Thus,  the 
procedure  of  using  the  Rayleigh  curves  with  a  median  signal-to-noise  ratio 
would  be  a  good  one.  On  the  other  hand,  if  one  is  interested  in  an  upper 
bound  on  detection  probability,  the  tail  of  the  density  function  must  be  con¬ 
sidered.  The  Rayleigh  curves  are  then  completely  unsu* table. 


Table  1.  Detection  Threshold  (Bias)  Level,  Y 

r. 


Detection  Threshold  (Bias)  Level,  Y, 
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Fig.  2,  Nonfluctuating  Target  (Marcum),  p  =  i  and  n  =  I0b 
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Fig.  3.  Nonfluctuating  Target  (Matcum)i  0-1  and  n  =  10 
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Fig.  4.  Nonfluctuating  Target  (Marcum),  p  =  1  and  n-  10^ 
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Fig.  5.  Case  1,  p  =  2,  n  =  Hr,  with  N  as  a  Parameter 
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Fig.  50.  Case  2,  n&  10  ,  N  =  1000,  with  p  as  a  Parameter 
Versus  Median  Signal-to-Noise  Ratio 
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Fig.  51.  Case  1,  p  =  2,  n=  10°,  Comparison  to 
Swerling  Curves 
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Fig.  52.  Case  1,  p  =  16*  n  =  10  ,  Comparison  to 
Swerling  Curves 
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Fig.  53.  Cate  2,  p  =  2,  n=  10  ,  Comparison  to  Swerling 
Curves  and  Normal  Approximation 
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Curves  and  Normal  Approximation 


-68- 


AVERAGE  S/N,  dB 

24  18 15 12  9  6  3  0  -3  -6 


lllll  illllllllllldllllllllllhlllllllPllllllllllllllllllkllUIIIIIIIIII 

lifliniiiaiiiiHiitUHiiimii  taaiiifluaiiiiiaifliiiiiiiuuifiiHiiiiii 

I  nan iii aaar aiaiiaaiaii  ii  !■■■■•*• ■■■■■■■■•■  '■■•■■■■  iit-iaiiiiii 


I  aaSiiiSa  laaaaaaaaacaiaaaa  ■■■■■•■>  liuiiiuNiiiBiikiaiiiatiiiiMinai 
I  aaaaa  aaa.  a  ■•■■■■•mi  aaaaa  ■■■■■  aat  a&  ■■■■■a  aaaaaaaaaa  9k  'iaaaaa  1  aa-»aa  saa&i 


CASE  2 

p  =  8 

n s  10® 


111131 
naail 

»iii 

■aMiitaii  1 

■■■■•■  25*5  | 

■aaaaa^agaaa  ■■■■•■■•  ■  “  ■■■■' 22!£i!535I?f? -2555S. 1 25222  SS5 52  SSSSSSSSSSSS SSSir  152  I 

[■■■■■•■a ■■■■■■■■■■■■■.'■  ■■■■*•■■■■  *«■■■■■■  '■■■■■!  ■«*■*■■■■■ ■■■■■■■aaa  ■■■■■■■■■■ 
IMMMIaaaiaaik  ■■■aaaanaiaat  naMaiaaiaiaaaaiiaaiaaaaa  maaiMHaganaiaiiiaiNiiaaal 
laaaaaaaiii ■■■•■*■*■■  iiaaaaaaaauia  taaiia i«aaa aaaaa aaaa»  •■■■•  ia**aaaaaa ■•«■■■■*■■  ■■■•*■■■■■  ■ 
■aaaaaaaaaa  )■•■  aaaSag a  Ja**« j*««a  ia *•■  aaay  4 ja«iaa«M laaaa  waaa  aaaajaaa*  aaaaaaaaaa  aaaaaaaaaa  ■ 

iiiiiijliiriiiiiiiitiiiiiiiiiiliiiiiiiaiiiiijiiiiiiii;iiiiiaiiiiiilHii|jiiiiiiiijiii 

hiiiiiiiiiu.iimiiitiiiinmiiiiniiiuiiiiifiiiiiiiiii[iiiii:iiiiiiiiiiiiiiimiiiiiiil 

iMIllllllllMHIUlMIIIIIIIIIHMIL'aillfilllHUIIHIIIIviniHlinilllllllllllllllllHl 
lik’iiiiiiiiiUkiiiiM iinmm '■riliiiiiiiiiiiiilimiiiml 
liiiMaiaHiiniHiaiiiraiiMH  jniiHi  '■■■UMiHMnniiii*aliiii.iiMiniiMiHaMi»iHH| 
lia^BaaaaBaiasai'aaaaaaaraiBaaaaiaaaaiagkaagaaBaaaaaaaaBaaaaaaaraa i^Baaaaaaaaaaaaaaaaaaaaaaal 

I  aikaaiaiai  ■  *■**  aaaaa  aaiai aaaaa  naaiiaii  liiaaaiaaniaaaaaaMiaiit  «•>  vaaaaaaaaaaaaaaaBaaaiaaaal 
liaa  i*  iaaaa  aaaaaaaaaa aaaaaaaaa.  aaaaa aaaa* aaaaaaaaaa  *aua«  aaaaa  aaaaaaaaaa  aaaaaaaaaa  ■ 

laaaLaaiaaiaaaaai  aaaa  aaa-ja  aaaaa  iaaaa aaaa*  ***** aaaaa  aaaaa aaaa a  aaaaa  aaaaa  aaaaa  aaaaa  aaaaa  aaaaa  I 

Jbaia  aHHi  a . .  iBMiiaaaa  <**•<■  aaaaa  aaaaaaaaaa •••  'a  u««a  aaaEaaaaaaaaaaiaaasafl 

■•aa  i  aaaaa  aaaaa  a  taai  aaacaaaaaa  aaaaa  aaaaa  a  at.  a  aaaaa  ■••MiiMi  a**t  a .  i«ii  aaaaaaaaaa  ■*•#■**•**■ 

^■•1  m  m «  m m mmm m  ■mmm  mmmm  mmmmmm  angaa  mmmmm  i  ■  Bi ■  m  mmmmm  Maaaaa  aaBaa  ■*■**  a  ■’■■a  aiaa*a  •**■■■  aaaaaaaaaa  ■ 

la  a  a  a  a aaaaa  aaaaa a >  aaa  aaaa  a aaaaaaaaaa aaaaa  »  mmmm  mmmmm  aaaaaaaaaa  aaaa  -  «  »aaa  aaaaa  aaaaa  aaaaaaaaaa ■ 
aaaaa  aaaaa  a aCaaaavaa  aaaaa  aaaaa  •_  TJS^SaaaS  5553S  SS5  Ja  aaaaa  5^252  jjSaag  aa  jyjJaaaaa  a  jgaa* 

I  iaaaa nmHBnmiliuiM ■■■■■ ■■•»?■■■■■■*”■■■■■«« ■■■■■■■■■a  ■■■■■*»  !aaiiaiaia**ii**ii«*a*«l 

I  •■■■  '■  u«aa iaaaa  iiNikim iai. ti iaaaa  ■■■■■■aaaa  ■«■■•«. \  ■■■■■■■  ■■■■■  ■■■■■■■>■■ 1 

1 1  iiMiMiiiaanilaiHim^utiaaaaiii  'tiui  iijaiiiiaiaiiaaaiiaaa*ii\aaiiiiiiaiaaaiiiiiaiiil 
la  •*aaa«r«a ■••■•■■■•  '■•■■■raaas ■■aaaaaa^ iiaat  a *<■■■■ ■■■■aaaaaa iaaaa ■■  <  maaaauaai iaaaa ■■■■■  I 
I  a,  mma  mmm  *a  ammmmmmmmi  aaaaaiiaaa  Haaaaaaai  ^aaa  ikaiMiaaaa iaaaa  laaaaiai »  a aaaaaaaaaa aaaaaaaaaa  I 

I  *i  •■(•■■JIIIMIIHII  ■■■■  >aa.  Mill  Hik  1  ■iHlMiai  Maiiiaiii  I 

bM|Maiiiiiiiiai*ii .  aaaaa  jaaa iaaaa  aaaaa  i*^5*  «4aa< aaaaaaaaaa aaaaa  aaa*  >  aaaaaaaaaa  I 

L  | 

■aag  a  a  n*M«amoa«*a»  *aaaJ 


■aaaaaa  •a'*!*  i-jaajiaaaaa 

■aaaaaa law %■ tuaaaiaaaaa _  __  . 

aaaaa  taa^o  im*  aaaaaaaaaa  aaaaa  aaaa  w  i 
-a  v ni ■■■■•uiat 


_  _ _ _ MflBBI  »■■■■■■■*  >•*»  > 

I  •«<■  •■■•«* i  aaaa  aaaaa  •*.  aaaa***  aaaaaaaaaa  ••■a**  '  ia«  ■aaalaaaaiaaaaaaaaaa  «« 

|  *  aaaaa  ■  aaaa  aaaaa  iat  uiawai  aaaaaaaaaa  aaaaa  •  ■  aaaaaaaaaa  aaaaaaaaaa  -v  'a 

I  aaaa  umm  iaaaaaafiaaaa  lafiaKii aaaaa aaaaa aaaaaa*.  ai aaaaaaiaaa aaaaaaaaaa r  a 

I  aaaa.  BBBIBBBBIBB,  ■  aaa  «*  iiaia  ■■Sib  aaaaa  iai  immirn  ««ai<  ■■  <M  Biaaa  w  » 

I  «••>«.»  «  •  mmm  mmmmm  mm  mm  i  mmmmm  mmmmm  mmmmm  mmmmm  mmm  -m  mmmmm  mmmmm  mmmmm  mmmmm  m 

I  aaaaa  aaaaaa  iaaaa*** aaaa.  i**w«ia**«a***a*aa*ai*a-  a aa**iiaa«« aaaaaaaaaa a^ 

I  aaaaa  t  ■■  mm  *iraa  aaaaa  uiu  «am  aaaaaaaaaa  aaaaa  aaiv  ■a%V4M«a>«a"a  aaaaa  aw. 

|  aaaaaa  «aaa  uvUMiiaiaaM  <M*i laaaailaaa  aian >iil  ta*fTVY*laa  aaaaaaaaaa  a«.> 

■ - -  -  - mmrnx  aaa  i  aaaaa  aaaaa  aaaaa  aaaa  4  r  I  fl  I  ■  laa  aaaaa  aaaaa  aa  . 

_ _ aaaa  ia»n  aaaaa  mmmmm  mmmmm  aaaa  «  •  a  t'l'iaB  aaaaa  aaaaa  « a  .  .  _ 

I  iiMit*  >aa  aaaa  1  aaaaa  aaaaa  B4  «ml  laaaiaaU*  «aaa«  lal*'1 1  »  fntfcla*  aaaaaaaaaa  a  ar  »  \  aaaaaaaaaa  a  aaaa  I 
I  i«a(M«t  Biiiia-iiiiaaaiiiiu  <»  ■■■aaiaaii  «■••« aaa*r  ■  traaaaaaa  aaaaaaaaaa  i»i4  aaaa  aaaaaaaaaa  I 
I aaaaa aaa  <aaawar  aaaaa aaaaa i»4 aaaaaaaaaa  aw  a  la aaaaa aaaaa aaaaa  aaa-u  4  a aaa aaaaaaaaaa I 
I  aaaa  a  aaa*  a  iaaaa  aaaaa  aaaaa  aaa  m  md  |  al  aaaaa  aaaaaaaaaa  ia  w  aaaaaa  aaaaaaaaaa  aaar  a  «  aaiiMaii«*n| 
I  aaaaa  aaaa  ■  aaaaa  jaaaa  aaaaa  aaaa  «  rl  1V|  aaaaa  aaaaaaaaaa  .ia»  a  •  aaaaa  aaaaa  aaaaa  aaaa  ■»  b  aaa  aaaaaaaaaa  1 
I  aaaaaaaaaa  iMaaiiauaaaatiaai.  r|Aai«aaaa  anaaaaaaa  k.aa  'k  taaaaaiaaa  aaaaa  aa»a  m  aa  aaaaa  aaaaa  I 
I  ■  aaaa  aaaaaaaaa  a^aai  aaaaaaaaaa  i  1»a  a  aaua  aaaaaaaaaa  aaaa  'i  ^  aaaa  aaaa*  aaaa  a  aaaat  aa  'aaaaaa  aaaaa  f 
I  «aa«aaiaia a  'aaa araaa aaaaaaaflaa k  :aiiiiutk5iaiiiiiiM««.  a  muiaiamiiiimnii  aaaaaaaaflia  I 
- - - - aaaa  a.  iaa  aaaaa  aaaaa  aaaai  ■r'aaa  u  aaa  iaaaa  aaaaa  a - -  - - 


SS1SSSS a 


I aaaaa  aaaaa  aaaaa  aaa.  a  raaaa  aa 

| aaaiijljaa iaaaa aaaa_ |iaa« §a 


I ••■■■•■•«• la 
Initlflilfilfl 


■■■■aaa 
■ - 


_ Mil 

1  nan  ai 


■PHiilliaiiinik' 
■i|Vlii  ■•■■■  If  ■■«  »fi|i 

iaiiv  Mail  ■■liiiairiifli 

flu  i«t^i  1 1  in  iiitt  din 


uiumiiiiiirliaMit  miiiii  laianu u  '<■111  i*k  n 
k^iiaiaiiiiiiiMiiaiiaN^iiiiiiaiiiiMiiiiriim- 
II  Ik'HlllllllltlcMlIlUi « .Bill  UIIIIIIIIIi‘11  <•■■■■  I 


I  mmmmS  Mm  mmm  mm  mMm  mmmmm  m  mm  mm  m  mii  mmmmm  mmmmm  mmmmm  mmmmm  a»Baa«a>**«B|  mi  mm  m  m  mm  Mm  m  mm  mm 

I  ««Bh*  ammm*  mmmmm  mmmmm  ••«#«  b*«b>  mmmmm  mmmmm  mmmmmmmmmrnmmmmmtimmmmmmmmmfmmmmmmmmmmmMmmmmMmmmmmmm 
I  aaaaaaaaaa  iia*<iai«*lMi«4<i*l|iH>aa«k«iti  iaaaa  baa  a  a  i  mtkmm  aaaaa  Uil  a  aaaaaa  «  iaaaa  aaaaa 

|  ■■aaiaaailMaaa aaiai  itlaiailiiilifeMaMiialHiaiiiiaiaiiAMi  maa naai  iiaii iia a naati aaaaa 

■  ■■>■«■  ■■■  •■■■■■■■■■■■■Mia  ■■■■■•■«  |i|iifin>  iMinifii  mMiiiii  iiniaiiti  i«iHifein 


••■■« ■■■•lllll 

viiiiliift Ilf  1 

iHBiiiiinii  I 

liiinaniili  I 

MIMHilllliJ 


— ■■■■«■ 

m  t  \  jt.'tti  iiiiiii  niiiiiiii  iimi  inn 

ilMiit  inti  |iin  niiigmi 

liiiBiiliiiftiiif  niiiiiiii 


RAYLEIGH 

L06-N0RMAL 


f lillslllslllilll II lllliili I 


0.2  0.4 


as  1.0 


1.8  2. 


R/Ro 


Fig.  55.  Case  2*  p  =  8,  ns  l(r.  Comparison  to  Swerling 
Curves  and  Normal  Approximation 
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Fig.  56.  Case  2,  p  a  16,  n  s  10  ,  Comparison  to  Swerling 
Curves  and  Normal  Approximation 
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Fig.  59.  Case  2,  p  =  2,  a  =  10  ,  Versus  Median 
Signal -to -Noise  Ratio 
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Fig.  60.  Case  2,  p  =  16,  ns  10°,  Versus  Median 
SignaUto-Noise  Ratio 
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Fig.  63.  Log-Normal  Probability  Density  Function, 
p  :  4,  x  as  a  Parameter 
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APPENDIX 


Computational  Procedures 


A.  CASE  1 


The  required  result  is,  from  Eqs.  (1)  and  (3), 


pd {*)=X  fN<y»x)dy 


=  J  dyj^  fN(y/x)  f(x,x)  dx 


Interchanging  the  order  of  integration  yields 


PD(x)  =Jdx  f(x,x)  J  fN(y/x)  dy 


The  inner  integral  is  seen  to  be  precisely  the  probability  of  detection  ob¬ 
tained  by  Marcum.  Denoting  this  by  P'(x)»  we  have 


PD(x)  *  J  p'(x)  dx 


The  integration  is  made  numerically  tractable,  taking  advantage  of  the  normal 
distribution  of  In  x,  by  means  of  the  substitution 


s  * 


In  x  -  In  x 

V*o 


(A-l) 
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T 


Then 


CD  2 

P  D(x)  =  /  e“Z  Px(z)dz 

**  mmCO 


(A-2) 


where  Pj(z)  is  P'(x)  under  the  transformation  of  Eq.  (A-l).  The  integration 
in  Eq.  (A-2)  is  readily  performed  by  means  of  a  numerical  quadrature  inte¬ 
gration  formula  based  on  the  Hermite  polynomials  [17].  As  in  Marcum  [5], 
P'(x)  is  computed  using  a  Gram-Charlier  series  with  the  Edgeworth  grouping 
of  terms  to  represent  fM(y/x).  Thus, 


fN(y/x)  ~7=  *l(t)  s  g(t) 

i=0  V^2 

where  t  =  (y  -  yJ/^TjT  Cli^  is  the  ^  central  moment  of  fN(y/x)], 


04(t)  =  e'*  U  H.(t) 

V^T  1 


and  H^t)  is  the  Hermite  polynomial  of  degree  i.  The  coefficients  Cj  are 
given  by 


(A-3) 


1=  -IT-/.H, 


(t)  g(t)  dt 


a  V-JL'  H4y5f)£N(y/x)  d* 


Substituting 


ti/2D  m  ,  j 


2mm!  (i  -  2m)! 


m=  0 


yields 


[i/2]  i+m 


2mm!  (i  -  2m)! 


m=0 

[i/2] 


E 


m=  0 


^i-2m 

,m _ *  / •  ^ _ i ■ 

2  m!  (i  -  2m)!  M2 


The  M}  are  expressed  in  terms  of  the  moments  about  the  origin  v.  by 


(A-4) 


n=  0 


The  moments  are,  in  turn,  obtained  using  the  characteristic  function  as  a 
moment  generating  function.  Thus, 


From  the  results 
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) 

(A- 5) 

9  p*o 

5:.. 


vv'*>=ferl,/2''y'*w2v^> 


(A-6) 


and 


C(p)  = 


-Nx 


Nx/(p+l) 

C 


Marcum  obtains 

vi  =  (1il -“TTf^  1F1 <**• N*  -Nx>  (A-?> 

where  jFj  is  the  confluent  hypergeometric  function.  The  coefficients  c^  are 
computed  from  Eqs.  (A-7),  (A-5),  and  (A-4)  using  the  recursion  relationship 

1F1(-i,N,Nx)=  Wk-Rvr:2i:--  ^(-i+l.N.Nx) 

'UTT^I  lF1(-i  +  2.N,Nx) 


Finally,  P'(x)  is  obtained  as 

P'(x)*/"fN(y/x)dy 

Yb 
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where 


The  function  fN(y/x)  is  given  by  Eq.  (A-3),  so  that 


The  result  of  the  integration  is 

P'(x)  =  -  [c3*2(T)]  -  [c4*3(T)+  c6«5(T)] 

-  [c5*4(T)+  c?$6(T)+  c9f7(T)] . 

The  brackets  indicate  the  Edgeworth  grouping  of  terms  (cQ  *  1  and 
Cj  s  c2  »  0).  Following  Marcum,  all  of  the  curves  except  those  for  N  =  1 
we:  e  computed  using  only  the  first  three  Edgeworth  groupings  to  approxi¬ 
mate  P'(x). 

The  accuracy  of  the  calculation  is  limited  by  the  number  of  ordinates 
used  in  the  quadrature  integration  formula.  A  total  of  20  were  used,  i.  e. , 

pl(a)d»-.^w(»i)p1(»i)  (A-8) 

i*  1 

where  the  are  seros  of  the  20^  degree  Hermite  polynomial  and  the  w(*j) 
are  the  quadrature  formula  weighting  fectors  [17].  The  effect  of  using  a 
small  number  of  ordinates  is  to  introduce  error  in  the  points  on  the  proba¬ 
bility  of  detection  curves.  In  general,  with  20  ordinates,  the  computed  points 
generate  a  smooth  curve,  although  some  spread  is  noticeable  for  the  larger 
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values  of  p.  The  manner  in  which  the  spread  decreases  is  shown  in  Fig.  A-l 
which  gives  points  for  N  =  1,  n  =  10*°,  and  p=  16  obtained  using  10,  16,  and 
20  ordinates.  It  appears  from  the  manner  in  which  the  points  converge  as  the 
number  of  ordinates  increases  that  the  maximum  error  in  probability  is  of 
the  order  of  0.  01  and  occurs  for  low  detection  probabilities.  The  curves  of 
P'(x)  obtained  for  N  =  1  and  p  =  1  are  quite  smooth  and  agree  remarkably  well 
with  Marcum's  curve  except  at  very  low  probabilities.  The  curves  of  P^(x) 
for  Ns  1  were  obtained  as  discussed  under  Case  2.  Cases  1  and  2  are  the 
same  for  N  =  1. 

B.  CASE  2 

The  probability  density  function  of  N  integrated  signal  plus  noise  variates 
at  the  detector  output  where  the  fluctuations  are  independent  from  pulse  to 
pulse  is  given  by  the  N-foid  self-convolution  of  f  ^(y,x).  For  numerical  compu¬ 
tation,  the  Hermite  quadrature  integration  formula  in  Eq.  (A-8)  is  a  convenient 
method  to  calculate  fj(y,x).  But,  for  this  case,  the  integrand  is  equal  to 
fjty/x)/^  under  the  transformation  of  Eq.  (A-l).  The  function  f  j(y,x)  was 
computed  for  621  points  over  the  interval  0(0.  06)37. 2*  for  each  combination 
of  p  and  x  (values  of  x  were  chosen  at  3-dB  increments).  The  function  i^{y,x) 
is  obtained  by  convolving  f  j(y,  x)  with  itself.  At  this  point,  the  increment  in 
y  was  doubled  so  that  the  311  points  on  the  interval  0(0. 12)37.  2  correspond  to 
every  other  point  on  the  previous  interval.  The  function  over  the  interval 
37. 32(0. 12)74.4  was  assumed  to  be  aero.  Similarly,  f^(y,x)  was  obtained 

*The  interval  (0, 37.2)  was  divided  into  increments  of  0. 06. 


from  f2(y,  x),  fg(y,  x)  from  f^y.x),  and  f1Q(y,  x)  from  both  f^  (y,  x)  and  fg(y,  x). 
This  procedure  is  continued  to  obtain  f^(y,  x)  for  N  =  20,  40,  80,  100,  200, 
400,  800,  and  1000.  The  increment  in  y  was  doubled  (in  the  same  manner  as 
before)  after  10,  20,  100,  and  200.  The  probability  of  detection, 

pD(x)  =  1  -  JQ  b  fN(y.x)  dy 

was  computed  by  Simpson's  rule  for  numerical  integration.  Simpson's  rule 
was  also  used  in  the  convolution  integrals. 

The  question  of  accuracy  is  difficult  to  answer  because  of  the  difficulty 
in  obtaining  the  maximum  value  of  high  order  derivatives  of  f^(y/x)  and 
f^(y,  x)  either  analytically  or  numerically.  These  derivatives  are  needed  to 
fix  a  bound  on  the  numerical  integration  errors  [17].  However,  by  doubling 
he  initial  increment  in  y  and  using  only  half  as  many  points  for  fj(y,x),  the 
resulting  computed  probability  of  detection  seldom  changed  by  as  much  as 
0.  01  (the  greatest  differences  usually  occurred  at  low  detection  probabilities). 
About  1  hour  of  computing  time  on  the  SDS  920  computer  was  required  to  find 
flOoo(y.x).  At  first  glance,  it  would  seem  possible  to  shorten  the  computation 
time  by  an  approximation  method.  This  is  discussed  next. 

As  the  number  of  integrated  pulses,  N,  is  increased,  it  is  expected  that 
fN(y,  x)  will  approach  a  normal  density  function  with  mean  and  variance  equal 
to  N  times  the  mean  and  variance  of  f^(y,x)>  It  is  also  possible  to  include  a 
first-order  correction  term  [18]  as 
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N 


? » x ) 


(A-9) 


where  ?  =  (y  -  |iN)/aN- 

2 

The  quantities  ^tN>  ,  and  are  the  first  moment  and  second  and 
third  central  moments  of  fN(y,  x),  respectively.  The  probability  of  detection, 
provided  the  approximation  is  valid,  is  given  by 


(A- 10) 


where  S  =  (Yfe  -  uN)/aN- 

The  general  moments  of  fj(y,  x)  can  be  calculated  in  a  straightforward 
manner.  The  ktk  moment  about  the  origin  of  fj(y/x)  as  a  function  of  x  is 
given  by 


/*  ® 

Vk(x)=/o  ykfj(y/x)dy 

yk  e"^X+y^  IQ(2^xy)  dy 
s  k!  Lk(-x) 

where  L>k(x)  is  the  Laguerre  polynomial  of  the  kth  degree.  Hence,  the  kth 
moment  about  the  origin  of  fj(y,x)  is  given  by 

vk*jT  vk(x)  f(x,  x)  dx 


which  can  be  expressed  in  terms  of  the  moments  of  f(x,x).  Thus,  the  first 
few  moments  are  given  by 


n  n(n-l)—  n 

x  =  p  x 

Hence,  the  moments  of  fj(y,x)  can  be  expressed  in  terms  of  p  and  x.  The 
central  moments  can  be  computed  by  means  of  Eq.  (A-4).  Finally, 

UN  =  N[1  +  x] 

0^  *  N[1  +  2x  +  (p2  -  l)x2] 
m3N  =  N^2  +  +  Mp2  *  l£t2  +  (2  +  p^  -  3p2 )5c 3  3 


For  large  p,  the  coefficient  of  the  correction  term  in  Eq.  (A- 10)  is  given 
approximately  by 


A  necessary,  but  not  sufficient,  condition  for  the  convergence  of  Eq.  (A- 10) 
to  the  true  probability  of  detection  is  that  Eq.  (A-ll)  should  be  much  less  than 
unity.  Based  on  this  result,  N  »  16  for  p  =  4,  N  »  1000  for  p  =  8,  and 
N  »  75,  000  for  p  =  16. 

The  probability  of  detection  based  on  the  normal  approximation  (one 

8 

term  only)  is  shown  in  Figs.  47  -  50  for  p  =  2,  4,  8,  and  16  and  n  =  10  .  For 
p  =  2,  the  normal  approximation  is  accurate  for  N  =  100  and  1000.  For  p  =  4, 
it  is  accurate  for  only  N  =  1000,  and,  for  p  >  4,  it  is  not  a  suitable  represen¬ 
tation  for  the  actual  probability  of  detection.  The  first-order  correction  term 
is,  in  most  cases,  either  negligible  or  unrealistic  (i.e. ,  probabilities  larger 
than  unity). 
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